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ABSTRACT 

A precise and accurate determination of the Hubble constant based on Cepheid variables requires 
proper characterization of many sources of systematic error. One of these is stellar blending, which 
biases the measured fluxes of Cepheids and the resulting distance estimates. We study the blending of 
149 Cepheid variables in M33 by matching archival Hubble Space Telescope data with images obtained 
at the WIYN 3.5-m telescope, which differ by a factor of 10 in angular resolution. 

We find that 55 ± 4% of the Cepheids have no detectable nearby companions that could bias the 
WIYN V-band photometry, while the fraction of Cepheids affected below the 10% level is 73 ± 4%. 
The corresponding values for the / band are 60 ±4% and 72 ±4%, respectively. We find no statistically 
significant difference in blending statistics as a function of period or surface brightness. Additionally, 
we report all the detected companions within 2" of the Cepheids (equivalent to 9 pc at the distance 
of M33) which may be used to derive empirical blending corrections for Cepheids at larger distances. 

Subject headings: Cepheid Variables; Galaxies: Individual (M33) 



1. INTRODUCTION 

An accurate and precise measurement of the Hub- 
ble constant at the few-percent level imposes signif- 
icant constraints on the equation of state of dark 
energy and other cq smologically relevant parameters 
(jKomatsu et al.l 1201 If ). The next generation of sur- 
veys aimed at improving our understanding of dark en- 
ergy will benefit from an even tighter constraint on Hq 
(Weinber g et alll2012f) than the present bounds of 3.4% 
jRiess et al.ll201lD . 

Cosmologi cal applications of the Ex tragalactic Dis- 
tance Scale (jFreedman fc Madord 120101 ) primarily rely 
on the Period-Luminosity rel ation of Cepheid variables 
(hereafter the "Leavitt law" , iLeavitt fe Pickering! 119121) 
as the pr imary distanc e indicator. The upcoming Gaia 
mission (jPrustil 1201 1|) is expected to deliver a sub- 
percent calibration of t he Leavitt law in the Milky Way 
(|Windmark et al. 2011), which could in turn enable a 1% 
measurement of Hq if all sources of systematic error are 
properly accounted for. 

One of these sources of systematic error occurs when 
two or more neighboring (but not necessarily physically 
associated) stars fall within the same resolution element 
of an instrument and cannot be fit with separate point- 
spread functions (PSFs). This effect is commonly re- 
ferred to as blending and it is different from crowding 
or confusion noise, which results in improper PSF fit- 
ting and/or inaccurate background subtraction due to a 
very high stellar density. An extreme example of blend- 
ing in the absence of crowding is a Cepheid in a binary 
system located in a low-surface brightness environment. 
Blending will bias the measured flux of a Cepheid to- 
wards larger values, shifting the Leavitt law to brighter 
magnitudes and leading to systematically shorter dis- 

1 Current address: Gemini Observatory, Northern Operations 
Center, Hilo, HI, 96720. jchavez@gemini.edu 

2 Current address: Department of Physics, Mount Allison Uni- 
versity, Sackville, NB E4L 1E6, Canada. 



tances and larger values of Ho. Extreme blends can 
be readily identified by their effects on Cepheid colors 
and/or a mplitude ratios and such tests are routinely car- 
ried out (IPellerin fc Macnl 120111: IScowcroft etal] [20091: 
iMacri et al.l 120061 ) . However, low- level blends are un- 
likely to be identified by such cuts and may affect stud- 
ies of the metallicity dependence of the Leavitt law (an- 
other source of systematic uncertainty) since they could 
mimic the photometric changes expected from differences 
in chemical abundances. 

The Local Group galaxy M33 is a good tcstbed 
for studies of Cepheid systematics tha nks to its rel- 
ative proximity (D = 895 — 965 kpc, Bo nanos et al.l 
I2006t IPellerin fc Macrj 120111). m oderate inclination an- 
gle (i = 55°. IHo et all Tl997l ) and recent episodes 
of star formation which have resulted in large num- 
bers of Cepheids throughou t i ts disk (lHartman et all 
[20061: IPellerin fc Macrl l201~ll) . IScowcroft et all ([20090 
used M33 Cepheids to study the "metallicity effect" 
of the Leavitt law, motivated by the la rge abundance 
gradient inferred from H II regions (Zaritsk v et all 
119941: iMagrini et al 1 120071 120 10ft. However other stud- 
ies (lUrbaneia et al.l [20051 iBresolin et all 120101: iBresolinl 
1201 1|) have determined a much shallower abundance gra- 
dient, which would make the metallicity effect consider- 
ably harder to measure. 

The disk of M33 has been extensively imaged by the 
Hubble Space Telescope (HST) using the Wide-Field and 
Planetary Camera 2 (WFPC2) and the Advanced Cam- 
era for Surveys (ACS). The angular resolution of HST at 
optical wavelengths is 10-15 times better than the see- 
ing at a good site on the surface of the Earth. Thus, a 
comparison of HST and ground-based images of the same 
Cepheids in M33 can yield useful insights into the nature 
of blending for more distant galaxies observed only with 
Hubble. 

Previous studies of the influence of blends on the 
Cepheid Distance Scale, based on comparisons between 
ground-based and HST images of nearby galaxies were 
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TABLE 1 

HST OBSERVATIONS OF M33 USED IN THIS STUDY 



Field 


RA 


Dec 


Camera 


Filters 


Exp. time (s) 


Prop. 


N 


Comments 


Name 


(J2000, deg) 






it^ 


TT 1 




# 


Ceph. 




b4w 


OQ flQQ A 
ZO.0O04 


on 70 -1 o 

oU. / olo 


WFPC2 


555 


814 


400 


400 


5237 


3 


16 in P07 


d2a 


ZO.4U0O 


Qfl C1 1 1 

oU.olll 


ACS 


606 


814 


10414 


20828 


9873 


7 


U49 in S06 


d4w 


OQ A f^l O 

zo.4D iy 


Qfl 70*77 
OU. li) 1 1 


WFPC2 


555 


814 


4800 


5200 


5914 


5 


U49 in K02 


fla 


oq ns 


Qfl 

ou.ozy i 


ACS 


606 


814 


10414 


20828 


9873 


6 


M9 in S06 


flw 


OQ KQQfi 


OU.Ooyo 


WFPC2 


555 


814 


4800 


5200 


5914 


1 


2 in C01 


glw 


oq f^ri/io 


oU.0o4o 


WFPC2 


555 


814 


4800 


5200 


5914 


3 


R14 in K02 


g2w 


OQ Qc: 
Zo.oLoD 


Qfl «A A P. 
OU.0440 


WFPC2 


555 


814 


4800 


5200 


5914 


13 


10 in C01 


g6w 


OQ /1fi7n 
Z0.4D ( U 


OU.DD i 1 


WFPC2 


555 


814 


1600 


1200 


5464 


2 




gbw 


OQ /IVC^Q 


Qfl ftQ7Q 


WFPC2 


555 


814 


2240 


2120 


6640 


2 


WW1 in C99 


gdw 


oq f^ono 


Qfl A O 

ou.oo4ty 


WFPC2 


606 


814 


340 


400 


8059 


2 




gew 


OQ /lOKQ 


OU.DZOU 


WFPC2 


606 


814 


1600 


360 


8059 


5 


also prop. #8805 


hlw 


oq /inQn 


ou. / uyo 


WFPC2 


606 


814 


8800 


17600 


9873 


4 




h4a 


OQ /I 1 7fi 


Qfl RA QO 


ACS 


606 


814 


2160 


2160 


10190 


6 


Fl in S09 


h5a 


OQ Qfi/1/1 


ou.DOoy 


ACS 


606 


814 


2400 


2500 


10190 


7 


F2 in S09 


jlw 


OQ 

Z0.0004 


oU. 001 ( 


WFPC2 


606 


814 


1600 


360 


8059 


2 


also prop. #8805 


j2a 


OQ A 


ou.oy t o 


ACS 


606 


814 


6380 


6624 


10190 


8 


innermost in W09; Dl in S09 


j2w 


OQ f^QflC 


oU. 0004 


WFPC2 


606 


814 


300 


700 


8059 


2 




j4a 


23.5010 


30.5487 


ACS 


606 


814 


2400 


2500 


10190 


5 


F3 in S09 


j5w 


23.4733 


30.6000 


WFPC2 


555 


814 


2240 


2120 


6640 


2 


15 in C01 


jaw 


23.4735 


30.5632 


WFPC2 


555 


814 


2240 


2080 


6640 


2 


17 in C01 


jbw 


23.4523 


30.5314 


WFPC2 


555 


814 


2240 


2080 


6640 


4 




jew 


23.4697 


30.5199 


WFPC2 


555 


814 


2240 


2080 


6640 


1 




11a 


23.5497 


30.4560 


ACS 


606 


814 


2160 


2160 


10190 


1 


F5 in S09 


mlw 


23.3728 


30.6106 


WFPC2 


555 


814 


520 


460 


6431 


10 




m3w 


23.3535 


30.5758 


WFPC2 


606 


814 


1600 


360 


8059 


10 


also prop. #8805 


m4w 


23.3179 


30.5566 


WFPC2 


606 


814 


240 


400 


8059 


4 




m5a 


23.4161 


30.4834 


ACS 


606 


814 


21260 


26420 


10190 


4 


D2 in S09, 26 in C01 


m6a 


23.3592 


30.5109 


ACS 


606 


814 


2160 


2160 


10190 


15 


F4 in S09 


n3a 


23.2727 


30.4538 


ACS 


606 


814 


2400 


2500 


10190 


4 


F6 in S09 


qla 


23.3348 


30.3706 


ACS 


606 


814 


21260 


26420 


10190 


6 


D3 in S09 


q2a 


23.3068 


30.3919 


ACS 


606 


814 


2160 


2160 


10190 


3 


3rd outermost in W09; F7 in S09 


q2w 


23.3645 


30.3742 


WFPC2 


555 


814 


4800 


5200 


5914 


1 


9 in C01; C38 in K02 



Note. — C99 : lUha ndar ct al. (19 93); C01: |C handar ct al. (2001); K02: Kim ct al. (2003); P07: IPark , 
ISaraiedini et alj 1 12009) ; S09: ISan Roman et al.l H2009I) ; W09: Willia ms et al.l 112009') 



"Lei {2007); S06: 



carried out by IMocheiska et al.l (|2000D in M31, by 
IMocheiska et al.l (|200lD in M33. and bv iBresolin et al ' 
( 20051) in NGC~300. In the case of M33. IMocheiska et al 



We based our 



( 20011) used HST/WFPC2 i mages and the Cep heid sam- 
ple of the DIRECT survey (jMacri et al.|[200lT) . During 
the intervening decade there have been numerous addi- 
tional HST observations of M33 using both WFPC2 and 
the Advanced Camera for Surve ys (ACS), which enable 
us to study more Cepheids than IMocheiska et al.l (|200lD 
and, in the case of ACS, with greater depth and finer 
pixel scale. Furthermore, we rely on a new synoptic sur- 
vey of M33 (jPellerin fc Macril 120111 ) carried out at the 
WIYN 3.5-m telescope with more Cepheids and better 
a ngular resolution than th e DIRECT catalog. 

iPellerin fc Macril (|201ll ) carried out extensive simula- 
tions based on the M33 ACS images to quantify the pho- 
tometric bias due to crowding in their ground-based pho- 
tometry. Considering the range of magnitudes and sur- 
face brightnesses spanned by the M33 Cepheid sample, 
they found that crowding bias increased as a function of 
magnitude but did not exhibit a dependence on surface 
brightness. Our paper complements their study by quan- 
tifying the photometric bias due to blending for Cepheids 
in M33. 

We describe in ij^jthe data used in this paper and the 
photometry we measured; describes the method used 
to quantify the level of blending; we discuss the results 
in 2] and compare them to previous work in fj5] Our 
concluding remarks and suggestions for future work can 
be found in 36l 



DATA AND ANALYSIS 

analysis on the Cepheid sam- 
jy the M33 Synoptic Stellar Survey 
nl l2011l ). We identify these variables in 
HST images and search for companions unresolved in 
the ground-based data. We calculate blending statistics 
based on these companions. 



pie published by the M33 
(jPellerin fc Macril 



2.1. Cepheid Sample 

Our analysis is based on the sample o f Cepheids listed 
in Table 3 of IPellerin fc Macril (|2011l ). The ground- 
based observations and analysis are described in detail 
in that publication, which we briefly sum marize here. 
Data from the DIRECT survey of M33 (jMacri et alj 
12001 were combined with new images obtained at the 
3.5-m Wisconsin- Indiana- Yale-NO AO (WIYN) telescope 
with the Mini-Mosaic (MiniMo) camera to detect 563 
Cepheids ranging in period from 2 to 110 days. The typ- 
ical FWHM of the WIYN images was 0"75, sampled at 
a plate scale of O'.'28/pix. The photometr y and astrom- 
etry w ere calibrated using the catalogs of iMassev et al.l 
([20061) . 

2.2. HST Data 

We queried the Hubble Legacy Archive (HLA) and 
the Mikulski Archive for Space Telescopes (MAST0 

3 The HLA and MAST are operated by the Space Telescope 
Science Institute (STScI). 
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for HST images of M33 obtained with either WFPC2 
or ACS which h a d ove rlap with the WIYN images of 
iPellerin fc Macrj ([201 1[ ). We selected observations with 
multiple exposures to allow for cosmic-ray removal. We 
also required a minimum of 100 s of total exposure time, 
to ensure a depth that would enable the detection of faint 
companions around the Cepheids. We further restricted 
our study to fields that were imaged in V (HST filters 
F555W or F606W) and I (HST filter F814W). 

The H ST fields contained 149 (~ 25%) of the Cepheids 
listed in IPellerin fc Macrj ([201 If ) . The locations of these 
fields are shown in Figure [Hand listed in Table [TJ The ta- 
ble also contains references to previously-published anal- 
yses of the data. Except for two ACS fields, all images 
were acquired on a single epoch and we therefore only 
have imaging of the Cepheids at a random phase within 
their pulsation cycle. 

The ACS images were reprocessed through the MAST 
On-The-Fly-Recalibration pipeline to apply the most up- 
to-date calibrations, while the WFPC2 images had al- 
ready been reprocessed using the final set of calibr ation 
frames in mid 2009 by STScI ([Gonzaga et al.ll2010h . We 
dow nloaded the reprocessed images and used MultiDriz- 
zle (jKoekemoer etal.1 12002;) to remove cosmic rays, cor- 
rect for geometric distortions in the cameras, and co-add 
multiple observations into master images. 

2.3. Photometry and Cepheid Search 

We performed point-spread functio n (PSF) pho tomc- 
try using DAOPHOT and ALLSTAR (|Stetsonlll987l ). We 
derived model PSFs using grids of a rtificial stars created 
with Tiny Tim ([Krist fc Hookll2004T) for the appropriate 
bandpasses, cameras and CCDs. We ran the FIND algo- 
rithm twice on each image, removing all stars found on 
the first iteration before proceeding to the second one. 
This increased the detection efficiency of faint stars, such 
as possible companions of a Cepheid. ALLSTAR was run 
one final time on the merged star list. Based on the ob- 
served luminosity functions, the photometry is complete 
to V ~ 25.5, 1 ~ 24.7 and V ~ 24.3, 7-23 mag for ACS 
and WFPC2, respectively. 

Instrumental magnitudes were converted to the HST 
VEGAMAG system using t he eq uations listed in Ap- 
pendix D of ISirianni et al.l (12005T) and th e coefficients 
listgdJnTablelO of iSirianni et all (|2005h and Table 2 
of lDolphinl (l2009h for ACS and WFPC2, respectively. 

Given the vastly different resolution and depth of the 
HST and WIYN images, the former had significantly 
larger stellar densities. Furthermore, the astrometric so- 
lution provided by the automat ed STScI pipeline is onl y 
accurate to a few arcseconds (jKoekemoer et all 120051 ) . 
We obtained a rough initial match between HST and 
WIYN images using the brightest few hundred stars in 
common. Once the gross astrometric offset had been 
removed, we matched the com plete star lists using DAO- 
MATCH and DAOM ASTER (|Stetsonlll~99l and refined 
the astrometric solution of the HST images. Cepheids 
we re then selected based on the coordinates tabulated 
bv IPellerin fc Macril ([20111 ). We visually inspected every 
Cepheid to ensure the star in the HST frame was indeed 
a match to the same star in WIYN image. This process 
helped to identify and correct a few erroneous matches 
where a faint star close to the Cepheid was originally 




Fig. 2. — Comparison of HST and WIYN images of three 
Cepheids in M33, illustrating different blending values. Left col- 
umn: WIYN V images; center column: HST V images; right col- 
umn: HST / images. Panels are 8" on a side and the black circles 
are 0.75" in diameter, equal to the WIYN FWHM. 

identified as the variable in the HST frame. Lastly, we 
estimated the disk surface brightness by averaging the 
background flux values reported by ALLSTAR for stars 
within 7" of each Cepheid. 

3. BLENDING CALCULATION 

We quan tify the level of b lendin g following the pre- 
scription of iMochejska et all (|2000| ). 

s F = Ecav/c (i) 

where Sf is the total flux contribution from the compan- 
ions relative to the Cepheid in filter F, fi is the flux of 
an individual companion star located within the critical 
radius and fc is the flux of the Cepheid. 

We calculated the values of S separately for V and 
I, using a critical radius of 0'.'375 which is the aver- 
age value of the half-width at half-maximum (HWHM) 
of the WIYN PSF. We only include companions that 
contribute 4% or more of the flux of a Cepheid in or- 
der to provide a conservative estimate of the blending 
value. This cu t-off was adopted by [ Mochcisk a et al.1 
(2000), although lMocheiska et all ([200J) raised it to 6%. 
In practice, stars with fi ~ 0.05/c (or Amag~ 3.25) 
are near the completeness limit of the ACS images rela- 
tive to the faintest, shortest-period Cepheids, which have 
V ~ 22.5,7 ~ 21.5 mag. 

In the case of Cepheids present in both ACS and 
WFPC2 images, we calculated blending values using the 
ACS data given its finer spatial sampling and increased 
depth. In the case of Cepheids present in multiple fields 
obtained with the same camera, we gave preference to 
the image with the deepest exposure time. If the ex- 
posures were of similar depth, we averaged the Cepheid 
magnitudes and the S values. 
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V-I (mag) 

Fig. 3. — Color-magnitude diagram of M33 Cepheids (in blue) 
and companions within 0'.'375 (in red), contributing more (filled) 
or less (open) than 4% of the Cepheid flux. Black dots are used to 
plot 3.5% of the stars detected in the ACS frames with I < 26 mag. 

Figure[3]shows a comparison of HST and WIYN images 
for three Cepheids with different values of S. Each panel 
is 8" on a side, centered on a Cepheid. Circles with 
radii of 0'.'375 (typical WIYN HWHM) are drawn around 
the variables. The Cepheids were chosen to show the 
range of blending values, from Sf = (top row) to Sf ~ 
0.6 (bottom row). The left column shows the WIYN V 
images, while the center and right columns show the V 
and I HST images. 

The photometry and blending values are listed in Ta- 
ble [2] For each Cepheid , we list the ID and period from 
iPellerin fc Macrl (pOll . the V magnitude and its un- 
certainty, the value of Sv and its uncertainty, and the 
corresponding information for the I band. Additionally, 
we tabulate the V and I surface brightness values and 
the designations of the WIYN and HST fields where each 
Cepheid is located. The uncertainties in Sf values are 
calculated by propagating the reported ALLSTAR pho- 
tometric uncertainties through Eqn. [TJ HST field codes 
are based on the field name listed in the first column of 
Table [TJ followed by a letter to identify the camera ('a' 
for ACS, V for WFPC2). 

Figure [3] shows a color-magnitude diagram of the 
Cepheids and all companions located within the critical 
radius. As a reference, we also plot 3.5% of all the stars 
with / < 26 mag detected in the V and I ACS frames. 
The companions span a broad range of colors and magni- 
tudes, but most are associated with the red giant branch 
and the red clump. These findings are not directly appli- 
cable to all Cepheid hosts, since different star-formation 
histories will alter the relative contributions of the upper 
main-sequence and the red giant branch. 

We used the HST star lists obtained in §2.31 to tabu- 
late all companions within a 2" radius of each Cepheid, 
presented in Table |3l Companions are labeled using the 
Cepheid ID from Tableland are numbered in increasing 



order of radial distance from the variable. We list the x-, 
y- and radial distance from the Cepheid (in arcseconds), 
the V magnitude and uncertainty, and the I magnitude 
and uncertainty. Some companions were only detected 
in one band. 

This extended dataset can be used for a variety of fu- 
ture studies. For example, comparisons of HST data 
with ground-based observations of M33 at different an- 
gular resolutions can be easily carried out by selecting 
the appropriate critical radius. Likewise, the sensitivity 
of blending values to the faint-companion cutoff limit can 
be explored. Lastly, suitable scaling of fluxes and angular 
separations can yield simulated HST images of Cepheids 
in similar environments out to D ~ 35 Mpc, at which 
point 2" at the distance of M33 would be equivalent to 
the angular resolution of HST in the V band. 

4. RESULTS 

We find mean blending values of Sv = 0.096 ± 0.015 
and Si = 0.083 ± 0.013 and median values of zero for 
both bands. Figure [4] shows cumulative distributions of 
blending values, while Figures [5] and [6] show the distri- 
bution of blending values as a function of period and 
surface brightness, respectively. Table |4] lists the frac- 
tions of Cepheids which meet several blending criteria as 
a function of period and surface brightness. We calcu- 
lated the uncertainty in each fraction using the binomial 
distribution approximation, 

<f) = y/f(l ~ 1W (2) 

where / is the fraction value and N is the number of 
Cepheids meeting a particular set of criteria. We cross- 
checked the validity of this approximation by performing 
100,000 bootstrap resamplings with replacement, which 
yielded the same uncertainties. Lastly, we tested the 
sensitivity to outliers in the distributions by performing 
the same number of jacknife resamplings, keeping 90% 
of the original sample. The derived fractions remained 
stable at the 2% level. 

The fraction of Cepheids with no blending is 
marginally lower (~ la) for Cepheids with P < 10 d than 
for ones with P > 10 d. Such a trend might be expected 
because the shorter-period, less luminous Cepheids can 
be affected (at a fixed flux ratio) by a larger fraction 
of disk stars. However, the difference vanishes when 
comparing the statistics of Cepheids affected at the 10% 
level. There is no significant difference in the statistics of 
Cepheids located in areas with "high" or "low" surface 
brightness. 

We also examined the effect of blends on the color of 
the Cepheids by calculating the value of Sy — Si for all 
Cepheids with non-zero values of either Sv or Si. The 
resulting histogram, presented in Figure [Jj shows that 
most blends do not appreciably change the color of the 
Cepheids: (Sv - Si) = 0.03 ± 0.27. 

5. COMPARISON WITH PREVIOUS WORK 

IMocheiska etaLl (|200ll ) analyzed WFPC2 images of 
M33 Cepheids discovered by the DIRECT project. We 
recalculated our blending values using the parameters 
adopted in that paper: a critical radius of 0'/75 and 
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TABLE 2 
Cepheid properties (Abridged) 



ID 


P 


V <T V 


Sv 




I VI 


Si 


<t Si SB V SB 7 


Field 




(d) 


(mag) 




(mag) 




(mag/d") 


HST WIYN 



J013332 


,36+302819. 


,8 


2 


.689 


21 


.635 





.038 


0. 


000 





.000 


21 


.008 


0. 


032 





.000 


0. 


.000 


21 


.92 


21 


26 


111 5 db 


mOf 


J013324 


20+302248. 


9 


2 


.695 


21 


.620 





.038 


0. 


.073 





.005 


21 


.369 





029 


0. 


.155 


0. 


.008 


21 


.51 


20 


.96 


qla 


mOm 


J013316 


.88+302157. 


.9 


3 


.187 


22 


.737 





.106 


0. 


316 





.087 


21 


.839 


0. 


055 





.000 





.000 


21 


.55 


21. 


.01 


qla 


mOm 


J013309 


.08+302354. 


.6 


3 


260 


22 


.215 





.067 


0. 


.000 





.000 


21 


.750 


0. 


043 





.000 


0. 


.000 


21 


.94 


21 


.39 


q2a 


mOm 


J013329 


.48+303614. 


3 


3 


760 


21 


.520 





.102 


0. 


000 





.000 


20 


.893 


0. 


125 


0. 


.000 


0. 


.000 


21 


.41 


20 


.56 


m 1 w 


mOj 


J013312 


.30+302355 


7 


4, 


.053 


21 


.457 





.072 


0. 


089 


0. 


.041 


20 


.787 


0. 


062 





.000 





.000 


21 


.87 


21 


.33 


q2a 


mOm 


J013322 


.10+303731. 


.9 


4 


.584 


20 


.768 


o 


.087 


0. 


000 


0. 


.000 


20 


.520 


0. 


082 





.000 


0. 


.000 


21 


.72 


21 


.03 


h.5a 


mOj 


J013428 


.61+304820 


7 


1 


.634 


22 


.243 


o 


.059 


0. 


000 


0. 


.000 


21 


.389 


0. 


085 





.000 





.000 


21 


.89 


20 


.45 


b4w 


m03 


J013350 


.94+303117. 


,4 


4, 


.767 


21 


.559 


o 


.060 


0. 


.158 


0. 


.015 


20 


.927 


0. 


035 





.564 


0. 


.025 


21 


.33 


20 


46 


i Viw 

J uvv 


mOc 


J013320 


.60+303458. 


,4 


4 


889 


21 


090 


o 


.137 


0. 


.123 


0. 


.024 


21 


.004 


0. 


128 





.000 


0. 


.000 


21 


.53 


20 


.85 


rn3w 


mOk 


J013319 


.10+303305 


6 


4, 


933 


21 


.445 





.048 


0. 


000 


0. 


.000 


20 


.931 


0. 


068 





.000 





.000 


21 


.30 


20 


.78 


ni4w 


mOk 


J013307 


.46+302512 


.6 


5 


.166 


21 


.949 





.037 


0. 


000 


0. 


.000 


21 


.250 


0. 


033 





000 





.000 


21 


.90 


21 


.37 


q2a 


mOm 


J013328 


.77+303753. 


,4 


5 


.176 


21 


.196 





.054 


0. 


163 


0. 


.016 


20 


.954 


0. 


065 





.364 


0. 


.050 


21 


.29 


20 


,40 


m 1 w 


mOj 


J013342 


.17+304841 





5 


177 


21 


.321 





.102 


0. 


040 


0. 


.005 


20 


.961 


0. 


.047 





.063 





.003 


21 


.76 


21 


.14 


d2a 


mOb 


J013341 


.39+304736 


,4 


5 


190 


21 


.837 





.051 


0. 


148 





.013 


21 


.152 


0. 


049 





.152 





.009 


21 


.69 


21 


.05 


d2a 


mOb 


J013316 


.08+302051. 


,8 


5 


204 


21 


.642 





.065 


0. 


000 





.000 


21 


.055 


0. 


058 





.000 


0. 


.000 


21 


.48 


20 


.91 


qla 


mOn 


J013405 


.98+303454. 


.0 


5 


.281 


21 


.034 





.048 


0. 


047 





Oil 


20 


.539 


0. 


024 





.000 





.000 


21 


.02 


20. 


.38 


i2w 


mOc 


J013321 


.78+303234. 


9 


5 


.282 


21 


.685 


o 


.298 





398 


0. 


.119 


20 


.858 


0. 


038 


0. 


.000 


0. 


.000 


21 


31 


20 


78 


ni4w 


mOk 


J013308 


.37+302607. 





5 


.314 


21 


.774 





.057 





000 


0. 


.000 


21 


.151 


0. 


055 


0, 


.000 


0. 


.000 


22 


.07 


21 


.39 


n3a 


mOl 


J013408 


.64+303754 


7 


5 


.316 


21 


.287 





.058 


0. 


000 


0. 


.000 


20 


.693 


0. 


089 





.110 





.021 


21 


.08 


20 


.20 


g2w 


mOd 


J013325 


.71+302138. 


.7 


5. 


.571 


21 


.538 





.045 


0. 


.076 


0. 


.007 


21 


.272 


0. 


029 





.128 


0. 


.010 


21 


.54 


21 


.01 


qla 


mOm 


J013321 


.51+302132. 





5 


.683 


21 


.555 





.066 


0. 


074 





.015 


20 


.974 


0. 


037 





.092 


0. 


.018 


21 


.54 


21 


.00 


qla 


mOm 


J013332 


.14+303002. 


3 


5 


.791 


21 


.803 





.108 


0. 


207 


0. 


.015 


20 


.981 


0. 


030 





.067 





.004 


24. 


.53 


23 


.07 


tii 6a 


mOf 


J013429 


.09+303735. 


.1 


5 


.828 


21 


.782 





.034 


0. 


.157 


0. 


.007 


21 


.584 





074 





.170 


0. 


.014 


21 


.81 


21 


.14 


fla 


m05 


J013425 


.11+303541 


.4 


5 


.868 


20 


.809 


o 


.033 


0. 


000 





.000 


20 


.667 


0. 


081 





.000 





.000 


21 


.80 


21 


.12 


fla 


m05 


J013350 


99+303156 


4 


5 


.889 


21 


,403 


o 


.091 


0. 


254 





.061 


20 


.762 


0. 


113 





.417 





.146 


21 


.14 


20 


.25 


irw 


mOc 


J013407 


,94+303831. 


,2 


5. 


.902 


21 


.060 





.028 


0. 


.000 


0. 


.000 


20 


.796 





072 





.144 


0. 


.024 


21 


11 


20 


.24 


gdw 


mOd 


J013359 


.89+303800 





5 


.997 


21 


.204 





.057 


0. 


000 





.000 


20 


.698 


0. 


065 





.000 





.000 


20. 


.52 


19. 


.83 


gew 


mOd 


J013349 


.55+304743 


,7 


6 


010 


21 


.718 





.031 


0. 


000 


0. 


.000 


21 


.161 


0. 


.078 





.000 





.000 


21 


.39 


20 


.53 


d4w 


mOb 


J013333 


.08+304230. 


.8 


6 


015 


21 


.802 





.065 


0. 


.000 


0. 


.000 


21 


.165 


0. 


053 


0. 


.000 


0. 


.000 


21 


10 


20 


.55 


hlw 


mOc 


J013350 


.67+303445 


.9 


6 


031 


21 


.063 





.064 


0. 


235 





.016 


21 


.013 


0. 


034 





.612 





.021 


21 


.20 


20 


.51 


j2a 


mOe 


J013325 


.62+303510 


.8 


6 


111 


21 


.010 





.040 


0. 


000 





.000 


20 


.659 


0. 


.084 





.000 





.000 


21 


.41 


20 


.71 


m3w 


mOk 


J013328 


.37+303730. 


.9 


6 


.116 


21 


.238 





.057 


0. 


.000 


0. 


.000 


20 


.621 





069 


0. 


.000 


0. 


.000 


21 


.30 


20 


.39 


mlw 


mOj 


J013353 


.70+303152. 


.0 


6 


.276 


21 


.709 





.103 


0. 


491 


0. 


.048 


20 


.788 


0. 


056 





.041 





.003 


21 


.53 


20 


.92 


j4a 


mOe 


J013346 


.45+304740. 


.5 


6 


333 


21 


.272 





.065 


0. 


.100 


0. 


.015 


20 


.557 


0. 


085 


0. 


.000 


0. 


.000 


21 


.31 


20 


,15 


d4w 


mOb 


J013401 


.87+304148. 


3 


6 


333 


21 


.714 





.035 


0. 


.000 


0. 


.000 


21 


.128 





071 





.369 


0. 


.033 


21 


.07 


20 


.16 


glw 


mOc 


J013400 


.58+303630. 


.6 


6 


.346 


20 


.846 





.021 


0. 


000 





.000 


20 


.335 


0. 


048 





.000 





.000 


20 


.68 


20 


.06 


gew 


mOd 


J013351 


,82+303310. 


.9 


6 


.504 


22 


.100 





.093 





000 


0. 


.000 


21 


.279 


0. 


041 





.417 


0. 


.175 


21 


.19 


20. 


35 


jaw 


mOe 


J013329 


.22+303136. 





6 


.530 


21 


.020 





.049 


0. 


.000 


0. 


.000 


20 


.496 


0. 


.042 


0. 


.000 


0. 


.000 


22 


,88 


22 


26 


m6a 


mOk 


J013315 


.71+303319. 


3 


6 


.557 


21 


,464 





.161 


0. 


084 


0. 


019 


20 


.615 


0. 


052 





.084 





.015 


21 


.45 


20 


.96 


m4w 


mOk 



TABLE 3 

Stars found within 2" of Cepheid variables (Abridged) 





ID 




Ax 


Ay 

(") 


Ar 


V 


TV 

(mag) 


I 




°7 


J013302. 


03+302553 


,4-cOOl 


-0.112 


0.010 


0.112 


24.021 


0. 


.434 










J013302. 


.03+302553 


4-c002 


-0.109 


0.221 


0.247 


25.333 





.134 


25 


.139 





.115 


J013302. 


.03+302553 


,4-c003 


0.333 


-0.132 


0.358 


26.490 





.091 










J013302. 


03+302553 


,4-c004 


-0.306 


0.260 


0.402 


26.442 





.190 


25. 


.226 





.338 


J013302. 


03+302553 


,4-c005 


0.375 


0.203 


0.427 


25.891 


0. 


126 


24. 


.051 





.085 


J013302. 


.03+302553 


,4-c006 


0.454 


-0.033 


0.455 


25.141 





.119 


24 


.446 





.099 


J013302. 


03+302553 


,4-c007 


0.047 


0.567 


0.569 


27.163 


0. 


113 










J013302. 


03+302553 


4-c008 


0.128 


0.652 


0.664 


27.007 


0. 


081 










J013302. 


.03+302553 


,4-c009 


0.109 


-0.743 


0.751 


26.728 


0. 


.069 


26 


302 





.096 


J013302. 


.03+302553 


.4-cOlO 


-0.790 


-0.128 


0.800 


22.802 





.071 


21 


.326 





.048 


J013302. 


.03+302553 


.4-cOll 


-0.787 


-0.223 


0.818 


25.424 


0. 


.247 










J013302. 


03+302553 


.4-c012 


0.069 


0.870 


0.872 


27.152 


0. 


145 










J013302. 


03+302553 


.4-c013 


0.877 


0.051 


0.878 


27.355 


0. 


143 










J013302. 


.03+302553 


.4-c014 


0.855 


-0.279 


0.899 


26.387 





133 


25. 


761 


0. 


.105 


J013302. 


.03+302553 


.4-c015 


-0.853 


0.323 


0.912 


27.702 





143 










J013302. 


03+302553 


4-c016 


0.628 


0.738 


0.969 


25.763 


0. 


101 


25. 


536 


0. 


.100 


J013302. 


.03+302553 


.4-c017 


0.289 


-0.941 


0.985 


26.810 





.114 


26. 


.077 





.119 


J013302. 


03+302553 


.4-c018 


-1.016 


0.074 


1.019 


24.672 


0. 


.065 


23. 


.880 


0. 


071 


J013302. 


.03+302553 


,4-c019 


0.588 


0.858 


1.040 


27.127 


0. 


113 










J013302. 


.03+302553 


4-c020 


-0.842 


0.663 


1.071 


26.431 





.096 


26. 


.194 





.111 


J013302. 


.03+302553 


,4-c021 


0.512 


-0.966 


1.093 


27.062 


0. 


111 


26. 


.481 


0. 


.094 
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Fig. 4. — Cumulative distributions of the blending values in the V and I bands (right and left panels, respectively). The solid lines 
represent the entire Cepheid sample while the dotted and dashed lines denote the short- and long-period (P < 10, P > 10 d) Cepheids, 
respectively. 



' 1* * * i 



t 



Period (days) 

Fig. 5. — Blending values as a function of the period of the Cepheid. 
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Fig. 6. — Blending values as a function of the sky background. 
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Fig. 7. — Distribution of the "color" 
to their Cepheid. 



of the companions relative 



TABLE 4 
Blending statistics 



Blending 


Sub- 


N 


Blending 


criteria 


level 


sample 




CMP12 


M01 








% 


% 


S v = 


all 


149 


55 ±4 


30 ±4 


S v = 


P < lOd 


72 


50 ±6 


22 ±5 


S v = 


P > lOd 


77 


57 ±6 


35 ±5 


S v = 


Sv < 21.4 


71 


61 ±6 


30 ±5 


S v = 


T, v > 21.4 


78 


47 ±6 


28 ±5 


S v < 0.1 


all 


149 


73 ±4 


45 ±4 


S v < 0.1 


P < lOd 


72 


74 ±5 


43 ±6 


S v < 0.1 


P > lOd 


77 


70 ±5 


46 ±6 


S v < 0.1 


S v < 21.4 


71 


76 ±5 


44 ±6 


S v < 0.1 


Y.V > 21.4 


78 


68 ±5 


45 ±6 


5/ = 


all 


149 


60 ±4 


30 ±4 


Si =0 


P < lOd 


72 


56 ±6 


19 ±4 


Si = 


P > lOd 


77 


61 ±6 


38 ±6 


Sj = 


Sj < 20.7 


74 


60 ±6 


24 ±5 


Sj = 


Sj > 20.7 


75 


57 ±6 


33 ±5 


Si < 0.1 


all 


149 


72 ±4 


41 ±4 


5/ < 0.1 


P < lOd 


72 


71 ±5 


33 ±5 


Si < 0.1 


P > lOd 


77 


71 ±5 


47 ±5 


Sj < 0.1 


Sj < 20.7 


74 


66 ±5 


32 ±5 


Si < 0.1 


Sj > 20.7 


75 


76 ±5 


48 ±5 



Note. — CMP12: this work; M01: IMocheiska et al.l pOOlT ) 

a companion flux cutoff of 6%. The results are tabulated 
in the rightmost column of Table 4. We also compared 
the individual blending values measured for 33 variables 
in common in V and 28 in /. As seen in Figure |8l there 
is good agreement with (ASf) — —0.02 ± 0.13. 

The statistic s deri ved using the criteria of 
IMocheiska et al.l (|2001[) are in excellent agreement 
with the values presented in their paper. For example, 
the fraction of Cepheids with Sy < 0.1 becomes 45 ±4%, 
compared to their value of ~ 43 ± 5% (inferred from 
their Fig. 4 and Table 2). We also obtain identical values 



Fig. 8. — Comparison of blending values for Cephe i ds in com- 
mon (found in WFPC2 images) with Mochciska ct al. (2001); our 
analysis was redone using their criteria. Filled (open) circles are 
used to plot the blending values in the V (I) filter. 

for the mean and median blending levels (24 ± 3% and 
13%) and reproduce the difference in blending statistics 
for "short" and "long" period Cepheids. Clearly, the 
differences between the two sets of values presented in 
Table 4 are due to the 2 x smaller critical radius adopted 
in our study, and emphasize the importance of angular 
resolution. 

We also compared the disk surface brightness values 
we d erived with those determined by (jMocheiska et al.l 
I2001f ) and found agreement at the level of 0.2 mag/d". 
We note that our surface brightness calculation is based 
on the average background le vel of the HST i mages 
within 7" of the Cepheid, while IMocheiska et al.l (|2001f) 
used the DIRECT ground-based images to calculate the 
sky in an annulus about 6" from the Cepheid. Regard- 
less of the method used to measure the background or 
the blending criteria adopted, there is little (if any) cor- 
relation between blending fraction and surface brightness 
f or the range of v alues c onsidered here. 

iBresolin et al.l (|2005) calculated blending statistics 
for a small sample of 16 Cepheids in NGC 300 using 
HST/ ACS images. They found a median value of 0% 
and an average value of 7%. Our results are consistent 
with their findings. 

6. CONCLUDING REMARKS 

We have presented a survey of Cepheids in M33 and 
their companions within 2", as resolved by HST with the 
ACS and WFPC2 cameras. We calculated the flux con- 
tribution of the companions when they are blended (un- 
resolved) in ground-based images with a seeing of 0'.'75. 
We find that more than half of the Cepheids in our sam- 
ple exhibit no blending at V and /, regardless of period or 
surface brightness. The majority of companion stars are 
located in the red giant branch and do not significantly 
alter the derived color of the Cepheids. 
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We plan to combine the ground-based photometry of 
iPellerin fc Macrj (|2011f) with the blending values derived 
in this paper to investigate possible biases in the determi- 
nation of distance moduli and "metallicity corrections" 
when using samples that lack such higher-resolution 
imaging. Additionally, our compilation of companions 
may be useful to derive empirical photometric bias cor- 
rections for Cepheids in more distant galaxies studied 
with the Hubble Space Telescope, provided the variables 
are located in similar environments to the M33 sample. 



JMC acknowledges support by the Department of 
Education through the GAANN Fellowship Program. 
AP and LMM acknowledge financial support through a 



Texas A&M University faculty startup fund. We thank 
Profs. Jianhua Huang and Lan Zhou for useful discus- 
sions on statistical techniques and the referee, Dr. Barry 
Madore, for his very helpful comments. 

Based on observations made with the NASA/ESA 
Hubble Space Telescope and obtained using the Mikul- 
ski Archive for Space Telescopes and the Hubble Legacy 
Archive (HLA) at STScI. STScI is operated by the As- 
sociation of Universities for Research in Astronomy, Inc. 
under NASA contract NAS 5-26555. The HLA is a col- 
laboration between STScI/NASA, the Space Telescope 
European Coordinating Facility (ST-ECF/ESA) and the 
Canadian Astronomy Data Centre (CADC/NRC/CSA). 

Facilities: HST (WFPC2, ACS), WIYN (MiniMo) 



REFERENCES 



Bonanos, A. Z., Stanek, K. Z., Kudritzki, R. P., Macri, L. M., 
Sasselov, D. D., Kaluzny, J., Stetson, P. B., Bersier, D., 
Bresolin, P., Matheson, T., Mochejska, B. J., Przybilla, N., 
Szentgyorgyi, A. H., Tonry, J., & Torres, G. 2006, ApJ, 652, 
313 

Bresolin, F. 2011, ApJ, 730, 129 

Bresolin, F., Pietrzynski, G., Gieren, W., & Kudritzki, R.-P. 

2005, ApJ, 634, 1020 
Bresolin, F., Stasiriska, G., VTlchez, J. M., Simon, J. D., & 

Rosolowsky, E. 2010, MNRAS, 404, 1679 
Chandar, R., Bianchi, L., & Ford, H. C. 1999, ApJS, 122, 431 
— . 2001, A&A, 366, 498 
Dolphin, A. E. 2009, PASP, 121, 655 

Freedman, W. L. & Madore, B. F. 2010, ARA&A, 48, 673 
Gonzaga, S., Biretta, J., & et al. 2010, HST WFPC2 Data 

Handbook (Baltimore: STScI) 
Hartman, J. D., Bersier, D., Stanek, K. Z., Beaulieu, J. -P., 

Kaluzny, J., Marquette, J.-B., Stetson, P. B., & 

Schwarzenberg-Czerny, A. 2006, MNRAS, 371, 1405 
Ho, L. C, Filippenko, A. V., & Sargent, W. L. W. 1997, ApJS, 

112, 315 

Kim, M., Kim, E., Lee, M. G, Sarajedini, A., & Geisler, D. 2002, 
AJ, 123, 244 

Koekemoer, A. M., Fruchter, A. S., Hook, R. N., & Hack, W. 
2002, in The 2002 HST Calibration Workshop : Hubble after 
the Installation of the ACS and the NICMOS Cooling System, 
ed. S. Arribas, A. Koekemoer, & B. Whitmore, 337 — h 

Koekemoer, A. M., McLean, B., McMaster, M., & Jenkner, H. 
2005, Demonstration of a Significant Improvement in the 
Astromctric Accuracy of HST Data, Tech. rep., Space 
Telescope Science Institute 

Komatsu, E., Smith, K. M., Dunkley, J., Bennett, C. L., Gold, 
B., Hinshaw, G, Jarosik, N., Larson, D., Nolta, M. R., Page, 
L.. Spergel, D. N., Halpern, M., Hill, R. S., Kogut, A., Limon, 
M., Meyer, S. S., Odegard, N., Tucker, G. S., Weiland, J. L., 
Wollack, E., & Wright, E. L. 2011, ApJS, 192, 18 

Krist, J. & Hook, R. 2004, The Tiny Tim User's Guide 
(Baltimore: STScI) 

Leavitt, H. S. & Pickering, E. C. 1912, Harvard College 
Observatory Circular, 173, 1 

Macri, L. M., Stanek, K. Z., Bersier, D., Greenhill, L. J., & Reid, 
M. J. 2006, ApJ, 652, 1133 

Macri, L. M., Stanek, K. Z., Sasselov, D. D., Krockenberger, M., 
& Kaluzny, J. 2001, AJ, 121, 870 



Magrini, L., Stanghellini, L., Corbelli, E., Galli, D., & Villaver, E. 

2010, A&A, 512, A63+ 

Magrini, L., Vflchez, J. M., Mampaso, A., Corradi, R. L. M., & 

Leisy, P. 2007, A&A, 470, 865 
Massey, P., Olsen, K. A. C, Hodge, P. W., Strong, S. B., Jacoby, 

G. H., Schlingman, W., & Smith, R. C. 2006, AJ, 131, 2478 
Mochejska, B. J., Macri, L. M., Sasselov, D. D., & Stanek, K. Z. 

2000, AJ, 120, 810 
— . 2001, ArXiv Astrophysics e-prints, astro-ph/0103440 
Park, W.-K. & Lee, M. G. 2007, AJ, 134, 2168 
Pellerin, A. & Macri, L. M. 2011, ApJS, 193, 26 
Prusti, T. 2011, in EAS Publications Series, Vol. 45, EAS 

Publications Series, 9—14 
Riess, A. G., Macri, L., Casertano, S., Lampeitl, H., Ferguson. 

H. C, Filippenko, A. V., Jha, S. W., Li, W., & Chornock, R. 

2011, ApJ, 730, 119 

San Roman, I., Sarajedini, A., Garnett, D. R., & Holtzman, J. A. 

2009, ApJ, 699, 839 
Sarajedini, A., Barker, M. K., Geisler, D., Harding, P., & 

Schommer, R. 2006, AJ, 132, 1361 
Scowcroft, V., Bersier, D., Mould, J. R., & Wood, P. R. 2009, 

MNRAS, 396, 1287 
Sirianni, M., Jee, M. J., Bem'tez, N., Blakeslee, J. P., Martel, 

A. R., Meurer, G, Clampin, M., De Marchi, G., Ford, H. C, 

Gilliland, R., Hartig, G. F., Illingworth, CD., Mack, J., & 

McCann, W. J. 2005, PASP, 117, 1049 
Stetson, P. B. 1987, PASP, 99, 191 

Stetson, P. B. 1993, in IAU Colloq. 136: Stellar Photometry - 

Current Techniques and Future Developments, ed. C. J. Butler 

& I. Elliott, 291 
Urbancja, M. A., Hcrrero, A., Kudritzki, R.-P., Najarro, F., 

Smartt, S. J., Puis, J., Lennon, D. J., & Corral, L. J. 2005, 

ApJ, 635, 311 

Weinberg, D. H., Mortonson, M. J., Eisenstein, D. J., Hirata, C, 
Riess, A. G, & Rozo, E. 2012, Physics Reports, in press 

Williams, B. F., Dalcanton, J. J., Dolphin, A. E., Holtzman, J., & 
Sarajedini, A. 2009, ApJ, 695, L15 

Windmark, F., Lindegren, L., & Hobbs, D. 2011, A&A, 530, A76 

Zaritsky, D., Kennicutt, Jr., R. C, & Huchra, J. P. 1994, ApJ, 
420, 87 



